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The form factors for the radiative transitions between charmonium mesons are investigated. We 
employ an anisotropic lattice using a Wilson gauge action, and domain-wall fermion action. We 
extrapolate the form factors to Q 2 = 0, corresponding to a real photon, using quark-model-inspired 
functions. Finally, comparison is made with photocouplings extracted from the measured radia- 
tive widths, where known. Our preliminary results find photocouplings commensurate with these 
experimentally extracted values. 
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1. Introduction 

The GlueX program at Jefferson Laboratory will produce hybrid mesons through photoproduc- 
tion, in contrast to the pion production mechanism generally employed in other searches. This pro- 
posal is supported by flux-tube model calculations, which show no suppression for conventional- 
hybrid photocouplings compared with those for conventional-conventional transitions[|l], |2[|3|]: 

a+y) ~ ^(100) keV, (1.1) 
p+y) ~ ^(1000) keV, (1.2) 
7T+7) = 230±60keV. (1.3) 

As a first examination of the validity of these calculations, we have embarked on a program to 
examine transition form factors in the charmonium sector. Here we expect both a flux-tube model 
and lattice calculations to provide reasonable descriptions of the physics. In this contribution, we 
focus on transitions between conventional charmonium mesons as a precursor to future studies 
of hybrid mesons. Many charmonium states below the DD threshold have measured radiative 
transtions into ligher charmonium states, and these transtions are reasonably well described in 
quark potential models [0]. 
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2. Formalism 

We write the transition matrix element between an initial state i and a final state / as 

(f{Pf,Xf)\V tl {Q,)\i(p u X i )) 

where pi and pf are the initial and final momenta respectively, A; and A/ label helicities, and is 
the electromagnetic vector current. For comparison both with experiment and with quark potential 
models, we express the transition form factors in a multipole expansion. For the transitions we 
consider in this paper, the expansion is as follows: 

2.1 T] c (0 ) form factor 

Whilst charge conjugation requires that this vanish, we can nevertheless examine the "form 
factor" by coupling the vector current only to the quark. This is the paradigm process, analogous 
to the pion form factor; continuum current conservation implies there is only a single form factor: 



(nc(Pf)\V' i (0)\r ]c (p i ))=f(Q 2 ) 



A 1 i M 



(2.1) 



where —Q 2 = q 2 = (Ef — E{) 2 — q 2 , with q = pf — p*j the three-momentum carried by the photon. 

2.2 J/y(1~) — ► 7T7c(0~) transition 

This is an M\ -transition, and likewise described by only a single form factor 

where X is the helicity label of the J/y. 



(r\c(p PS ) I V(0) | Y(Pv,V) = 2V(Q2) e* g %gsP?e y (Pv,A), (2.2) 
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J/Y 


XdO 


Xcl 


h c 


Lattice (MeV) 


2819(7) 


2917(7) 


3288(15) 


3401(29) 


3351(19) 


PDG (MeV) 


2980(1) 


3097 


3415 


3511 


3526 



Table 1: Determination of the spectrum, together with the corresponding PDG values [[$]]. 

2.3 Xco(0 + ) — ► y.//y/"(l~) transition 

This transition involves an electric dipole E\(Q 2 ) coupling to the transverse component of the 
photon, and a further form factor C\{Q 2 ) coupling to the longitudinal component: 

(Xco(Ps)\V^(0)\J/Y(py,X)) = 
Q.- l (Q 2 )[E l (q 2 ) [a(Q 2 )e»(p w ^)-s(p w ,l). Ps (PvPv-Ps- m2 vPs)} + 



1 m v s(pv,X).p s [pw-PsiPN+PsT -mlp^-m\p^\ (2.3) 
For real photons, only E\ is of interest, but in our calculation we will extract both form factors. 

3. Computational details 

The computations are performed in the quenched approximation to QCD on 300 configurations 
of a 12 x 48 lattice. We employ an anisotropic Wilson gauge action [Q], with a renormalized 
anisotropy ^ = a s /a t = 3. The temporal lattice spacing obtained from the static quark-antiquark 
potential is of 1 = 6.05(1) GeV. 

The quark propagators are computed using an anisotropic domain-wall fermion (DWF) action, 
with a domain-wall height M =1.1 and L5 = 16. We employ standard Gaussian smearing for the 
interpolating fields, using sources of various widths. For this exploratory study, we do not attempt 
a precise tuning of the charm quark mass to yield the physical m^. The charmonium spectrum 
obtained on our lattices is listed in Table [I]; in general, the whole spectrum is too light, reflecting 
the imprecision in our tuning. 

The three-point functions are computed using the usual sequential-source method, following 
the strategy outlined in ref. Sequential-source propagators are computed for both a final scalar 
Xco and a pseudoscalar r] c , placed at the mid-point of the lattice tt = 24, for both pt = (0,0,0) 
and pf = (1,0,0) where the momentum is expressed in the appropriate lattice units. The form 
factors are computed using the local, non-conserved vector current. The form factors are obtained 
from the three-point functions with the remaining amplitudes and masses extracted from fits to the 
two-point functions. The matching factor to the continuum, Zy, is determined by imposing charge 
conservation on the pseudoscalar form factor f(Q 2 = 0) = 1. We find a discrepancy of around 11% 
between the matching factor obtained for a pion at rest, and that for a pion at the lowest non-zero 
value of the lattice momentum p = (1,0,0), representing an important systematic uncertainty on 
our calculation. 
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Figure 1 : The green and red data points show the lattice determination of the r\ c "form factor" obtained using 
the spatial and temporal components of the electromagnetic current respectively. The band corresponds to a 
quark-model-inspired fit to the data, while the dashed line is the VMD expectation, as described in the text. 
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Figure 2: The left- and right-hand figures show plateaux in the determination of V(Q 2 ) at several Q 2 for 
Pf = (0,0,0) and pf = (1,0,0), respectively. 



4. Results 



The tj c "form factor" is shown in figure [lj for the case pf = (0,0,0). In order to describe the 
<2 2 -dependence of the data, we appeal to a non-relativistic quark model using harmonic-oscillator 



wave functions, characterized by a single parameter j3: 



f(Q 2 



(4.1) 



A single-parameter fit to the data yields j3 = 0.522(5) GeV, a value not uncharacteristic of that 
typically employed in quark-model charmonium wavef unctions; the fit is shown as the band in the 
figure. Also shown as the dashed line is the VMD form, using the lattice value of the J/y mass; 
VMD provides a notably poor description of the data, in contrast to the case of the pion[^]. 

Proceeding now to the transition form factors, we consider first the transition J/y — ► yr] c . 
The quality of the extraction of the form factor is illustrated in Figure || where we show the plateaux 
at various Q 2 using the spatial components of the vector current. 

The Q 2 dependence of the transition form factor is shown as the left-hand plot in Figure |3[ The 
systematic difference arising from the different determinations of Zy is manifest in the discrepancy 
between the values obtained with Pf = (0,0,0) and those with pf = (1,0,0); the right-hand plot 
shows the data with those for pf = (1,0,0) rescaled by this 11%, greatly reducing this discrepancy. 
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Figure 3: The left-hand plot shows the transition form factor J/y — > yrj c , for Pf= (0,0,0) and Pf = 
(1,0,0). The purple and yellow bursts are phenomenological values obtained from the PDG[^[], and using 
the transition rate to 00 respectively. The right-hand plot shows the same data with those corresponding to 
Pf ■ = ( 1 , 0, 0) scaled by the 1 1 % discrepancy in Zy ; the band is a quark-model-inspired fit to the lattice data, 
as described in the text. 

In order to extrapolate to Q 2 = 0, we once again employ a non-relativistic quark-model-inspired 
form 



v(Q z ) = v(oy 



(4.2) 



A combined fit to the p f = (0,0,0) and rescaled p f = (1,0,0) data yields j8 = 0.52(1) GeV and 
7(0) = 1.92(2). 

The radiative width is related to the photocoupling V(Q 2 = 0) by 



r(J/y -> i] c y) = a 



{m w + m nc ) 2 3 



|V(£ 2 = 0)| 2 , 



(4.3) 



whence we can obtain an "experimental" determination of the photocoupling using the measured 
radiative width. This is shown for the PDG value as the purple burst [|3|]; an alternative value 
obtained from the ratio of the product branching fraction of y — ► yr] c — ► y<p<p and the r\ c — ► <p<p 
branching fraction is shown in yellow. 

Finally, we turn to the XcO — ► yJ/W transtion. The two form factors E\{Q 2 ) and C\{Q 2 ) are 
shown as the left- and right-hand plots in Figure |J Using a quark-model form for the extrapolation 
to Q 2 = 0, we have for the electric-dipole form factor 



E 1 (Q 2 )=c\q(Q 2 )\e- Q2/ ^ 2 
where q is the three-momentum of the photon in the rest frame of the decaying XcO, 

_ (rn\ - m\) 2 + 2Q 2 {m 2 v +m 2 ) + Q 4 



4m 2 



(4.4) 



(4.5) 



with an analogous form for the longitudinal multipole C\{Q 2 ). The extrapolations are shown as 
the bands on the figure. Once again, we can use the measurement of the radiative width to obtain 
an "experimental" value for the physical multipole E\ (Q 2 = 0). This is shown as the purple burst 
using the PDG width|^|], and as the yellow burst using a recent CLEO determination^]. The 
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Figure 4: The left- and right-hand plots show the form factors E\(Q 2 ) and C\(Q 2 ) respectively for the 
radiative transition XcO — > YJ/W- The bands correspond to independent, quark-model-inspired fits to the 
lattice data at pf = (0,0,0) and pf = (1,0,0), as described in the text. The purple and yellow bursts are 
phenomenological values obtained from the PDG[|]|, and from a recent CLEO determination [^] 

corresponding values for the parameter j8 are 0.60(3) GeV and 0.47(1) GeV for the E\ and C\ 
multipoles respectively. 

In this work, we present a first study of radiative transitions in charmonium. While the results 
are preliminary, we find quark-model extrapolations yield values for the photocouplings in reason- 
able proximity to experimental expectations, with broadly consistent values for the wave-function 
parameter /3 . Computations of further form factors are in progress. A future goal is to extend this 
study to the light-quark sector relevant for JLab, and to transitions involving hybrid mesons. 

Acknowledgements 

This work was supported by DOE contract DE-AC05-84ER40150 under which the Southeast- 
ern Universities Research Association operates the Thomas Jefferson National Accelerator Facility. 
We are grateful for discussions with Ted Barnes, Frank Close, George Fleming and Jim Napolitano. 

References 

[1] F. E. Close and J. J. Dudek, Electroweak production of hybrid mesons in a flux-tube simulation of 



lattice qcd, Phys. Rev. Lett. 91 (2003) 142001 [hep-ph/0304243] 



[2] F. E. Close and J. J. Dudek, Hybrid meson production by electromagnetic and weak interactions in a 



flux-tube simulation of lattice qcd, Phys. Rev. D69 (2004) 034010 flhep-ph/ 0308 980 - 
[3] S. Eidelman et. al, Review of Particle Physics, Physics Letters B 592 (2004) 1+. 

[4] E. J. Eichten, K. Lane and C. Quigg, B meson gateways to missing charmonium levels, Phys. Rev. Lett. 



89 (2002) 162002 [htep-ph/0206018| 



[5] T. R. Klassen, The anisotropic wilson gauge action, Nucl. Phys. B533 (1998) 557-575 



[hep-lat/9803010] 



[6] Lattice Hadron Physics Collaboration, F. D. R. Bonnet, R. G. Edwards, G. T. Fleming, R. Lewis and 



D. G. Richards, Lattice computations of the pion form factor, hep-lat/041102 



[7] CLEO Collaboration, N. E. Adam, Branching fractions for psi(2s) to j/psi transitions, Phys. Rev. Lett. 



94 (2005) 232002 [tiep-ex/050302 



210/6 



